Collisional destructions of giant planets and rare types of meteorites 
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The probabihty of coUisions and destructions of giant planets at the various stages of planetary 
systems evolution is calculated. The flux of destructed planet fragments and the probability of their 
observations near the Earth are estimated. Of the particular interest is a case of the metastable 
metallic hydrogen fragments. The radio bursts, which can be generated by the collapsing magneto- 
spheres of the giant planets during mutual collisions, are also discussed. 
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I. INTRODUCTION 

In the recent years a lot of exoplanets — the planets 
outside the solar system in orbits around of the distant 
stars have been found. Most of the discovered planets 
belong to the class of giant planets, like Jupiter, because 
of the observational bias, but the solid-surface extrasolar 
planets of the Earth type are also detected. It is already 
clear that the planetary systems are not rare, but a typ- 
ical astronomical phenomenon. 

The planets in the solar and other planetary systems 
have being born in the gas-dust disks, resulting from the 
compression of protoplanetary clouds along the axes of 
common rotation. The terrestrial planets formed by the 
aggregation of dust particles and larger planetesimals. 
The giant planets formed, probably in the accretion of 
gas onto the solid seeds with the mass of a few Earth 
masses, or in the alternative model — from the gas con- 
densations due to the gravitational instability in the disk 
After the start of nuclear reactions in the central 
star, light elements remaining in the gaseous phase were 
carried out to the periphery by the stellar wind, with the 
exception of material which is already condensed into the 
protoplanets. 

At the various stages of planetary system evolution the 
planets can scatter gravitationally, collide, merge or de- 
struct, with the release of fragments of different sizes and 
compositions The catastrophic collisions are pos- 

sible especially at the stage of planet formation, when 
mutual scattering of planets are most often, and their 
orbits are rather chaotic and unstable. The process of 
planetesimal collisions and merges is considered as the 
mechanism of aggregation and growth of the terrestrial 
planets and the central stony cores of giant planets. The 
collisions with the following merges or destructions of the 
terrestrial planets were considered, for example, in [^-Q. 
Collisions and merges of giant planets may also be im- 
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portant factors in the early evolution of at least of some 
part of them. There are several mechanisms for the ap- 
pearance of overlapping planetary orbits and collisions: 
collisions at the early chaotic stage, convergence and in- 
tersection of the orbits in their slow change (migration), 
the resonant and three-body interactions. Finally, the 
convergence of orbits can occur in the decay of plane- 
tary system due to the mass loss of the hosting star Q, 
including the loss of massive shells or the central star su- 
pernova explosion. Immediately after the explosion, the 
planets fly apart from the star system center, and the 
inner planets can catch outer ones, having larger veloc- 
ity of the orbital motion. In this article we estimate the 
probabilities of these processes. 

The consequences of mutual giant planet collisions de- 
pend, in particular, on their internal structure, the ratio 
of planet masses, the relative velocity, the impact param- 
eter and angular momentums. For the terrestrial planets 
with solid surfaces, the effects of collisions were consid- 
ered in 6]. For the giant planets, the depth of interpen- 
etration of gas (or liquid) shells of two planets is very 
important, since the density in the outer layers increases 
rapidly inward. In the inelastic collision with a small rel- 
ative velocity, if the collision affected the large fraction of 
mass, the planets lose the kinetic energy of the relative 
motion and merge into the single body. We describe this 
process in some details below in the Section IIVI) . Other- 
wise, if the planets experienced a shallow tangential col- 
lision and then flew away, some of the planet substance 
can be released in a region of the compensated gravity. 
Collisions of giant planets must be accompanied also by 
the collapse of the planet magnetospheres, leading to the 
generation of the powerful bursts of radio waves available 
for registration. 

What fragments could remain after the partial de- 
struction of the giant planets in the collisions? External 
gaseous layers, clearly, quickly disperse. Heavy elements 
from the central cores of giant planets form a popula- 
tion of meteorites with properties similar to the ordinary 
meteorites, but possibly with some chemical anomalies. 
Special interest is in the fate of the intermediate regions 
of the giant planets, which arc probably composed of the 
liquid metallic hydrogen (MH) 10 12]. The possibility 
that hydrogen at high pressure can be in the solid and 
metallic phases was predicted in [l^, and the equation 
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of state of MH was discussed in [T^l ■ Review of the the- 
oretical works on the problem of MH is presented, for 
example, in [l^. Witnesses of the MH have been ob- 
tained in several experiments on the shock compression 
and the compression of hydrogen in the diamond anvil 
cells, but the results are so far controversial. Bibliog- 
raphy of the recent experimental work on the creation 
of MH can be found in the Introduction of the article 
[l6j . According to the theoretical calculations [Hi , 
MH can be metastable when the pressure is removed, 
but this conclusion is still under the debates (see details 
in the review [l^). The loss of stability is due to the 
sub-barrier quantum process of electron capture with the 
transition into the molecular phase, which occurs even at 
zero temperature. The quantitative characteristics of the 
metastable state were not clarified yet, for example, the 
lifetime of the sample depending on its size and environ- 
mental conditions is unknown. It is also unclear whether 
the MH transforms into the solid phase, when the pres- 
sure is removed and the sample is cooled from ~ 10'* K 
(the temperature inside the planet) to the temperature 
~ 10^ K (for orbit near the Jupiter in the solar system), 
corresponding to the balance of radiation heating and 
cooling. 

In view of potential practical significance of the pres- 
ence of MH in the space, we estimate the number of such 
fragments on the assumption of the metastability of MH, 
leaving a lifetime of the metastable state as a free pa- 
rameter. We will also discuss the characteristics of the 
passage of pieces of the MH in the Earth's atmosphere. 

The prehistory of the formation of the solar system 
is important for the probability of detection of the frag- 
ments (conventional rocks or hydrogen) in space from the 
destroyed planets-giant and the fragments collision with 
the Earth. Clash of the giant planets in the solar sys- 
tem at the early stage is the most promising source of 
the presence of MH fragments in the solar system. The 
Sun, according to the content of the metals, is the star 
of the second generation (not counting the pre-galactic 
stars) formed in the gas cloud left over from the explo- 
sion of a star of the previous generation. The formation 
of the solar system could began from the compression of 
the cloud by the shock wave from the near (at a distance 
of ~1 pc) supernova. Fragments of the MH could remain 
from the collision of planets belonging to the previous 
exploded star, they could also fly from the explosion of 
a nearby supernova, or from the more distant regions of 
the Galaxy. 



II. RATE OF GIANT PLANET COLLISIONS 

At the initial stage of the planetary system formation, 
the planets have rather chaotic orbits Q , they often over- 
lap, and the planets collide and merge, or are partially 
destroyed. In particular, the impact hypothesis explains 
the formation of Moon by the collision of a large body 
with the Earth. In addition, a possible explanation of the 



particular configuration of the planetary orbits in the so- 
lar system is obtained if one assume that there was a 
fifth giant planet, which has been ejected from the solar 
system in the scattering flQ*! . This may indicate that the 
chaotic stage with the mutual intersections of the orbits 
influenced not only the terrestrial planets but the giant 
planets too. It should be noted that the majority of the 
detected exoplanets are very close to the hosting stars 
than the giant planets in solar system, so the mechanisms 
of the formation of the different planetary systems may 
be very different from each other, at least quantitatively. 
While it is difficult to judge how often the chaotic orbits 
stage realizes, and whether it was in our solar system. 

In the number of papers devoted to the formation of 
planetary systems, the effect of "orbits migration" — the 
slow orbit change due to interaction (exchange of the an- 
gular momentum) with the gaseous disk — was consid- 
ered (see e.g. 0,[23l- For orbits with different ellipticity 
and different inclination to the ecliptic, the migration can 
lead to the intersection of the orbits and to the final col- 
lision of planets. A possible case of the rapid migration 
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two stages: at first, the planet app roaching of to the star, 
and then fiying away from it [2J]. In particular, Jupiter 
in its history first approached the Sun to a distance of 
about 1 AU, and then returned to its present position. 
The unstable motion of the planets may occur after the 
dispersing of the protoplanetary gas cloud, which has the 
stabilizing effect on the orbits. 

The cross section for two spherical massive bodies col- 
lision is ii 
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where ri^2 and mi_2 are the radii and masses of bod- 
ies, Wroi is their relative velocity. The first term in the 
brackets corresponds to the purely geometrical collision 
cross section and the second term takes into account the 
gravitational focusing. The unit in the parentheses, as 
a rule, can be neglected. For example, for planets with 
masses 7111^2 — M and radii ri_2 — R of the Jupiter and 
Dioi — 13.1 km s~^ (the orbital velocity of Jupiter), the 
second term in the brackets is equal to ~ 1 1 . If we con- 
sider the process of the destruction of the inner layers 
of a gas-giant planets, then the ri 2 in the formula ([T]) 
should be considered as the radii of these layers, but the 
mi 2 remain equal to the total masses of the planets. For 
example, calculations of the internal structure of Jupiter 
indicate that the MH occupies the region from 0.2 to 0.8 
of the radius of the planet, and the internal 20% of the 
radius fills the rocky core [24I. To estimate, we assume 
that the protoplanetary disk enclosed N giant planets in- 
side the characteristic radius L and thickness 2H. The 
thickness of the disk can be estimated from the charac- 
teristic values of the angle a of the orbit's inclination to 
the ecliptic (or to the median plane of the protoplane- 
tary disk) H ~ aL. For most of the planets in the solar 
system the characteristic value of this angle is ~ 1°, for 
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example, for Jupiter a = 1.03°, and for Saturn a = 2.5°. 

Consider first the case of chaotic orbits. We do not 
take into account the dependence of the quantities on 
the distance from the star and we restrict the estimations 
only by the average values. The characteristic time-scale 
of the planet's collision, having the concentration n2 ^ 
N/{nL^a), is 



for the case of the orbit migration 
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2.4 X 10^ yr. 
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where the last numerical estimate corresponds to the set 
of the values: a = 2.5°, L = 9 AU, i^roi — 13.1 km s~^, 
TV = 5, i? = i?,7 « 7 X 10* km, M = Mj w IQ-'^Mq- 
I. e., the parameters of the planets were normalized to the 
mass and radius of Jupiter, and the radius of the disk was 
putted equal to the radius of the Saturn orbit. It should 
be noted that during the evolution of the planetary sys- 
tem, its radius could vary significantly. For example, in 
the so-called "Nice model" of the solar system formation 
the initial configuration of the protoplanetary cloud was 
more compact than in the modern solar system [2]| . 

The motion of the planets is gradually becoming a reg- 
ular, so that the convergence of the Keplerian orbits at 
the later stage can not be considered as accidental. Reg- 
ularization of the orbits occurs, including the their cir- 
cularization due to interaction with the central star and 
the gaseous disk Q. If the planetary orbits were ellipti- 
cal initially with the possibility of intersections (like Nep- 
tune and Pluto), then later they can be circular, and the 
probability of collisions is reduced. Nevertheless, many 
of the observed extrasolar planets' orbit are not circular, 
but have large enough eccentricities instead, that may 
be due to inefficiency or the circularization and mutual 
scattering of the planets in their encounters. 

Consider now the case when the orbits are almost sta- 
bilized, but experience the migration with characteristic 
time-scale tj^ig ~ 10^ yr. The characteristic velocity of 
the migration is Vnug ^ L/{2t^ig) ~ 2 x 10^'' km s^^. 
We introduce the parameter of randomness 
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where Toib is the average orbital period of the planets, 
and Tp — ((7coii/7r)^/^ is the impact parameter of collision. 
The condition of randomness rj > 1 implies that during 
the orbital period the planet comes out from the volume 
swept by the cross section ([T|) , and the collisions are inde- 
pendent in each period. The time of the collision in this 
case is given by ([2]). Otherwise, if ?7 < 1, the planet at the 
next period will pass partially through the region of space 
where it was previously, and therefore the probability of 
the collisions is suppressed. If 77 < 1 the collision cross 
section sweeps the volume of space ~ tWmig2rp27rrorb dur- 
ing the time t, where Torb L/2 is the average radius of 
the orbit. Using, as above, the characteristic parameters 
of the protoplanetary disk, we estimate the collision time 
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Since the quantity Q is order of magnitude greater than 
the characteristic time-scale of the migration, the prob- 
ability of the collisions of the planets at the stage of the 
migration will be ~ 1/10. 

After finishing of the chaotic phase and stabilization 
of the orbits the collisions are possible due to the slow 
changes of the orbits under the influence of weak per- 
turbations. The probability of collisions during the sub- 
sequent slow evolution (after the dispersing of the gas 
disk) is small, and the calculations require the numerical 
simulations. 

On the basis of smallness of the characteristic times- 
scales in the simple estimates © and it is reason- 
able to assume that there were fs ^ 0.1 — 1 collisions 
of giant planets such as Jupiter for each planetary sys- 
tem during its formation. In the Galaxy, during the time 
^ 10^° yr, there were 10-^^/s destructions (by the num- 
ber of stars). In result, the rate of giant planet destruc- 
tions in the Galaxy is ^ 10/s yr~^. 



III. COLLISION AFTER SUPERNOVA 
EXPLOSION 

The closer the planet to the star, the greater its or- 
bital velocity. If the star exploded as a supernova, and 
thus completely or largely lost weight, then the planets 
fly away freely and there is a possibility that the inner 
planets will catch up the outer planets and collide. The 
explosion of the core-collapsed supernova in the case of 
low initial mass leads to the star total disintegration. 
More massive pre-supernovaes form the compact rem- 
nants. Thermonuclear supernovae of type la occur in the 
binary systems, and in this case the existence of extraso- 
lar planets around double stars is also possible, like the 
recently discovered extrasolar planets Kepler-16 b, -34 b, 
and -35 b. If the explosion is due to the flow of mate- 
rial from the companion star into the white dwarf, the 
companion remains after the explosion. If the alternative 
scenario is realized with the merge of two white dwarfs, 
these dwarfs are completely destroyed by the explosion. 
Loss of the planetary system is possible with the explo- 
sions of any of these types of supernovae, if the loss of the 
central mass will be more than half of its initial value, 
because in this case the initial circular velocity of the 
planets exceeds the escape velocity at the same distance 
from the star: GKh/r > 2GMf/r for Mf/Mi < 1/2. 

In the case of the slow mass-loss with the characteris- 
tic time much longer then the orbital period, the circular 
orbit of the planet will expand, but still remain circular 
due to the conservation of the adiabatic invariant. Oth- 
erwise, in the case of the sharp mass-loss, the circular 
orbits of the planets become elliptical with the possibil- 
ity of mutual intersections. Thus, any sharp loss of stars' 
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mass leads to the expansion of the orbits and increase of 
their elhpticity, which can cause colhsions. 

Planetary system after the supernova explosion, or af- 
ter the change of its structure, fly apart during the single 
passage through the system ~ 1 yr, so the probabil- 
ity of the collision can be estimated using the results of 
the Section ini as ~ ^i/^coir Taking into account that 
supernova explosions happen in our Galaxy about once 
every tgN = 20 — 50 years (we take 35 years for the es- 
timation), one can expect that the rate of the collisional 
destructions of the the giant planets is 

PsN - ii/(OsN) - 10-^ yr-\ (5) 

Over the lifetime of the Galaxy, ^ 10^° years, about ^ 
10'^ exoplanets was destroyed in this way, which is several 
orders of magnitude less than the number of destructions 
at the chaotic stages, that was discussed in the Sectionllll 

Increasing the radius of the orbit increases the like- 
lihood of the ejection of planets or other bodies from 
the system under the influence of tidal forces from of 
the neighbour stars in the rare encounters. Only the 
tidal mechanism contributes to the full release of the 
fragments from the star, because the total mass of star 
and shells dropped by the supernova remnant is approxi- 
mately equal to the initial mass of the star, so the pieces 
can not depart too far from the ejected shell, remaining 
gravitationally bounded to the system. Otherwise, due 
to the tidal forces, they can leave the system and migrate 
to the Galaxy. 

In collisions with relative velocity exceeding the 2nd 
escape velocity the complete destruction of the planets 
can occur with the large ejection of the matter. Such 
events can occur after the supernova explosion, if at least 
one of the planets was very close to its star and has the 
large orbital velocity. Only a rocky planet or a rocky 
core of the former giant planet could be close to the star. 
KIC 05807616 is the example of the system in which the 
giant extrasolar planets have lost the outer shells while 
increasing the radius of the star, and the only rocky cores 
orbit the star [2^ after the reverse compression. The shell 
ejected in the explosion of a supernova should also have 
a strong impact on the giant planets. The span of the 
shell and the subsequent emission of the supernova can 
cause destruction of the planet's outer layers under the 
mechanical and thermal effects. 



IV. MERGE AND FLY AWAY OF PLANETS AT 
COLLISIONS 

In order for the collisions of the planets to destroy 
them, the relative collision velocity (which is an order 
of magnitude of the orbital velocity Worb) must be larger 
than the escape velocity V2 on the surface of the planet 
(for the solid planets the excess is 2-3 times). For a planet 
with the parameters of Jupiter, orbiting a star with the 
parameters of the Sun, the orbital radius where this con- 
dition is satisfied, is r = GMq/vI =0.25 AU, because for 
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Figure 1; Relative mass loss l\M/M in dependence of the 
radius of penetration r (in units of planetary radius R) of 
the giant planets with the internal structure similar to the 
structure of Jupiter. The horizontal dotted lines indicate the 
AM/M, for which the collision touches the region of metallic 
hydrogen {AM/M ^ 1/80) or the merging of the planets 
occurs due to the loss of kinetic energy of their relative motion 
(AM/M ~ 1/20). In the intermediate region the fragments of 
metallic hydrogen can release. For comparison, dashed curve 
shows the mass loss for the model of Saturn. 



the Jupiter V2 = 60 km s~^. So, in the solar system at 
present there are no giant planets that can experience a 
collision with the complete destruction. The condition 
V2 < Vorh is valid for many exoplanets, which because of 
their proximity to the stars called "hot Jupiters". It is 
unlikely that in the solar system at an early stage of its 
evolution were the giant planets near the Sun. Therefore, 
this case is applicable probably only to other planetary 
systems. Due to the unfavourable temperature condi- 
tions near the star, a small piece of MH in such the envi- 
ronments is likely to quickly lose stability and evaporate. 

In the opposite case, V2 > Worb, the planets in the head 
on collision will merge and the most of the fragments 
will combine into the single large planet. Of the partic- 
ular interest are the tangential collisions of the planets 
when the planets fly away after the collision, having ex- 
perienced only a slight fracture of the outer layers. In 
the space between them in the region of the compen- 
sated gravitation part of the MH fragments may remain, 
which will be distributed over the different orbits around 
the star. Let the relative velocity of the planets before the 
collision is less than V2 velocity on the surface of planets. 
In the geometric intersection of the collision the matter 
stops, and its kinetic energy is partially converted into 
the heat, i. e. the inelastic collision of planets occurred. 
Part of the energy dissipated, but the total momentum 
and angular momentum with respect to an arbitrary cen- 
ter is conserved. The planets continue to fly apart, and 
the stopped substance is attracted by the gravity of the 
planets (due to the specified conditions of the velocities), 
roughly equally separated. Only in the region between 
the planets, where the attraction is compensated, part of 
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the substance remains. It loses its connection with the 
planets, and continues to move in the orbit that would 
have the center of mass of the pair of planets at the mo- 
ment of the collision. The released material, probably 
composed of many fragments of various sizes that fly 
close to, but at gradually diverging orbits. In all cases, 
the duration of the collision compared with the orbital 
period around the star can be neglected, assuming that 
the collision occurs at a single point of the orbit. 

Now we will find the condition for the release of MH 
from the interior of the planets when they collide. Cross- 
ing of the planets has the form of cylindrical surface that 
was cut in sphere that then collapses to the center un- 
der the influence of gravity. To estimate the captured 
material's mass we approximate the collision area by the 
plane cut of the sphere, whose lower flat area is close to 
the center of the planet at a minimum distance r. The 
mass AM of this figure is 



AM = -nR^ J dx{l - x^)p{Rx), 

r/R 



(6) 



where p{r) is the density of the planet in dependence on 
the distance till the center r, and R is the full radius of 
the planet. The density p(r) we take from the results of 
calculations of the internal structure of Jupiter, shown 
at Fig. 35 in the book fTTl. Fig. [T] shows the relative 
affected mass AM/M of the collisions in dependence on 
the minimum distances to the center. The area if MH 
will be affected if r/R < 0.8 (for the model of the internal 
structure of Jupiter). Thus, from Fig. [T] it follows that 
the minimum release of matter in this case is 



AM/M - 1/80. 



(7) 



The relative velocity v of the planets at the moment 
of the collision exceeds Urei due to the additional kinetic 
energy obtained during the approach: 
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However, when considering the fly of matter, it is more 
convenient to take into account the initial energy of the 
planets at large distances before the collision, when they 
still have a relative velocity Vrd- Let us assume for esti- 
mates that about half of the mass what have been lost 
during the collision of the planets then was pulled back 
by the gravitation of the planets. The energy loss on the 
attraction of the mass is [/ - G{AM/2)M/R. The plan- 
ets do not merge into a single body, if U does not exceed 
the initial kinetic energy of relative motion where 
the reduced mass fi = MiM2/{Mi^M2) = M/2, and the 
velocity Wrci ~ '^orb- For the collisions of the planets with 
the parameters of Jupiter, at the orbit with the radius 
''orb of Jupiter's orbit around the star with the mass of 
the Sun, this criterion gives 

AM/„,<_£_. 4X10- (A) (_-_)-,„) 



that is, AM/M J ~ 1/20. Thus, between © and © 
there is the mass interval in which the ejection of matter 
is possible. In this case the mass in the region of the com- 
pensated gravity can be estimated as Mh ~ Mj/lOO « 
1O~^M0. If this mass is homogeneously distributed over 
the volume V ~ of the protoplanetary disk, dis- 

cussed in section |TT1 the matter would have the density 
pi - Mh/V - 2.7 X lQ-14 g cni-3. 

The probability of the collisions with the release of MH 
can be estimated from the cross section ([1]), which, for the 
velocity condition under consideration, is proportional to 
the first power of (ri -t- r2). From (j9l) and Fig. [1] one ob- 
tains that the planets will merge with the probability 
68%, the release of MH will occur with the probability 
80 — 68 = 12%, and the remaining probability 20% cor- 
responds to the tangent collision which touches only the 
outer atmosphere without the direct influence onto the 
MH. 

In the region of the compensated gravity the liquid sub- 
stance is likely to have a very dispersed form, composed of 
fragments and droplets of various sizes. Thermal history 
and the fate of these droplets will be analyzed further. 
Fragments of the MH which were attracted by the planets 
will fall on the planet or will be revolve around the plan- 
ets as satellites. In the future, due to tidal interactions 
with other planets, these fragments will have chances to 
leave the planet and to go onto the independent orbits 
around the star. 



V. SWEEPING OF DEBRIS TO THE GALAXY 
BY TIDAL FORCES FROM NEARBY STARS 



Stars in the Galactic disk have peculiar velocities with 
the dispersion 20 — 50 km s""'^. Approaches and the grav- 
itational interactions of stars at distances ~ 1 — 10 AU, 
where their velocity vector is turned by ^ 7r/2, are very 
rare events occurring only once in 10^^ — 10^^ years, but 
such encounters destruct the planetary systems with the 
release of the planets, asteroids, and the other fragments 
into the interstellar space. Approaching at distances 
much larger than the size of planetary systems are more 
likely, but they only perturb the motion of the planets 
without causing an immediate release of the substance. 

We will speak about the solar system for definiteness. 
Let the star with mass 77i^, ~ Mq fly near the Sun with 
the impact parameter I — minimum distance in the case 
of a rectilinear trajectory. We assume that the star flies 
at a speed greater than the velocity at the outer orbits 
under consideration. Then one can use the impulse ap- 
proximation, which assumes that the perturbed body re- 
mained near stationary during the action of the disturb- 
ing forces. Note that the impulse approximation breaks 
down for the planets close to the Sun, because of the 
large characteristic frequencies of the orbital motion of 
inner planets, so the inner planets are protected by the 
conservation of the adiabatic invariant. We direct the 
axis z from the particle of matter to the trajectory of 



6 



the star at the point of the closest approach. For the 
straight-hne trajectory (impulse approxiraation) , the an- 
gle between the velocity of the star and the line con- 
necting the particle and the star is evolving as follows 
d(p/dt = — tiroicos^ (p/l. After the change of variable t to 
the (f) in the equation of Newton one finds 



Vrcll 
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(10) 



and the result of integration — extra velocity of the par- 
ticle in the direction z is given by 



2Gm^ 



(11) 



For the removal of matter from the solar system, not the 
absolute rate of the gain, but its increment with respect 
to the velocity of the Sun is important. Denoting by Vz 
the increase of the velocity of the Sun, one finds that the 
relative increment is 

I'z - =i = ^rcosV', (12) 

where r is the particle distance from the Sun, ip is the 
polar angle from the axis z in spherical coordinates. The 
increase in the velocity in the transverse direction has 
the same order of magnitude. We will not consider the 
exact configuration of the flight and orbit, and restrict 
ourselves only to the order of magnitude estimations, in 
which the total increase is 
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One must to distinguish the two cases: (i) the destruc- 
tion of the planetary system by the single close star pas- 
sage and (ii) the perturbation due to many distant tidal 
interactions with stars. In the first case, the condition 
Aui Worb gives l/r ^ V2(forb/t^roi)"^/^ ~ 1, i-e, as ex- 
pected, the passage should be close to the size of the 
planetary system in this case. The probability of such 
encounters during the time tg is 



Pi ^ Trr^n* Vicing ^ 2 x 10 



(14) 



where n* 0.12 pc^^ is the concentration of stars in the 
solar neighbourhood, Wrei ~ 30 km s~^ and tg ~ 10^° yr is 
the typical age of a planetary system (for the solar system 
tg ~ 5 X 10^ yr), and the radius of the orbit r ~ 10 AU 
was taken. As can be seen from the above estimates, 
the probability of destruction of the planetary system is 
very small, but in the whole galaxy ~ IQ-'^^Pi ~ 2 x 10^ 
such events could occur with the removal of the planets, 
as well as various pieces of matter into the interstellar 
space. 

We now consider the accumulation of many weak in- 
teractions due to the distant encounters. Since the incre- 
ments are not correlated, the increase of the energy has 



the diffusion form, the squares of the energy increases are 
summed. Using (jl3p we obtain 



d{Avf 

dr~ 



AG'^mlr'^n 



dvvf(v) f 2'Kldl 



(15) 



where f{v) is the velocity distribution of stars. In the in- 
tegral over the impact parameters I the close approaches 
give the main contribution, but the rarity of them is taken 
into account automatically in (jisp due to the smallness 
of n*. After integration we obtain the average growth 
rate of the planets velocity during the time tg-. 

Av - 2Gm^nJ'^tl/^/vlli^ - 0.1 km s^^ < Vorh- (16) 

Thus, in the typical cases, the tidal interactions do not 
destroy planetary systems, but only add the little chaotic 
velocity component to the planets. 

Tidal interactions, in particular, can be responsible for 
the large inclination « 17° of the orbit of Pluto to the 
ecliptic plane and its eccentricity. This feature of the 
orbit may be the result of one or more encounters with 
stars, under the assumption that Neptune and Uranus, 
during the interaction were on the other side of its orbit, 
and therefore did not experience the strong perturbation. 



VI. RADIO BURSTS 

The collisions of giant planets can generate the bursts 
of radio waves due to the destruction of their powerful 
magnetospheres such as Jupiter's magnetosphere. Bursts 
of energy can in many respects resemble the solar fiares 
with close characteristic time scales. 

Now we will estimate the energy fiux of the electro- 
magnetic burst. The magnetic field at the visible surface 
of Jupiter is about Bi ~ 10 G, and its magnetosphere 
extends for a distance of 100 Jupiter radii. The colli- 
sion excites the plasma turbulence and provides the gen- 
eral compression of the matter. For the one-dimensional 
compression, the magnetic field can increase till the value 
Bf ^ 10^(>^/10^)B^. The energy of the magnetic field 

^'^^'-^^^ 6xl0-(^)\rg. (17) 



E„ 



3 Stt V102 
This energy is released in a time At ^ Rj /v^ci ^ 
1.5 hours, and therefore the power is Pm ^ Em/ At ^ 
10'^^(>f/10^)^ erg s~^. If the burst occurred at the dis- 
tance Im, then the observed fiux is 



Prr 



10-1^ 1 



VloO 



10 kpc 



erg s ^ cm ^. 



(18) 

We assume that the emission of the energy in the 
plasma in magnetic field occurs at the cyclotron fre- 
quency uJc/2tt = eBf/{2Trmec) ^ 3 GHz, then the spec- 
tral density of the flow is 

Fm ^ ^ 3 X 10-^ {^] { TT-^ I Jy, (19) 



Vl02/ V 10 kpc 
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Figure 2: The rate of the radio bursts generated by the col- 
lisions of giant planets in the Galactic disk. The logarithm 
of the number of bursts per year is shown on the vertical 
axis with the observed flux on the Earth > Fm, measured in 
Jansky. Two curves show the two cases with the number of 
collision per each star fs =0.1 and fs = 1 during the time 
~ 10^" yr. 



that is enough to register it on the current radio tele- 
scopes, which have the typical sensitivity of ^ 10 ^Jy. 

Now we will find the distribution of radio bursts ob- 
served from the Earth, taking into account the struc- 
ture of the galactic disk. We assum that the collisions 
of the planets are uniformly distributed over the time 
tg - 10^° yr, and we use the usual exponential model of 
the surface density of the stellar disk 



(Tsir) 



-r/ro 



(20) 



VII. FRAGMENTATION OF PLANETARY 
MATERIAL AND THERMO-BALANCE OF 
DEBRIS 

The gravitation in the giant planets in large enough 
volumes plays a more important role than the cohesive 
forces in the liquid, surface tension etc. In [7] the study of 
the coUisional destructions of gravitationally-dominated 
bodies was performed in the model in which the bodies 
were composed by the densely packed spherical particles 
with gravitational attraction. Gravity gives some effec- 
tive plasticity to the substance of the planet, because the 
force of gravity tends to return the planet into a spherical 
state after the disturbance. Thus, in terms of strength 
of materials subject, a collision goes through the plastic 
deformation in the shear environment with the thermal 
energy then release — the transition of the kinetic energy 
into the heat, and the evaporation of some fraction of the 
mass. 

It can be assumed that by the time of the collision 
the areas of liquid MH in the protoplanets have already 
formed. MH which avoided the immediate evaporation 
will be distributed in the form of fragments (droplets) of 
various sizes. We take the spectrum of the debris in the 
form of the simple exponential function, and the index 7 
will be regarded as a free parameter: 



Nf{m)dm 



JWmax 



mm 



(22) 



The normalization is performed on the total mass that 
could release into the fragents, as it has been estimated 
in the Section ITVl 



M„ 



M„ 



mNf{m)dm = fhMh, 



(23) 



where Md = 8 x lO^^M© and tq — 4.5 kpc. In this equa- 
tion the distance of the star from the center of the Galaxy 
r(Z, 6) = [P + — 2lrQ cos 9]^^^ is expressed through its 
distance from the Earth I, where tq = 8.5 kpc is the dis- 
tance of the Sun from the galactic center, 9 is the angle 
between the directions to the star and to the center of the 
Galaxy. Then the rate of bursts with the signal greater 
than Fm is 



N{> Fm) = / d9 



dlfjst-^cJs{r{l,9))/m, 



(21) 



is shown at Fig. [21 where lm{Fm) is the inverse of 
and fs = 0.1—1 was evaluated in the SectionlTTl Thus, one 
sees that the radio bursts are available for the registration 
if they could be distinguished from the backgrounds or 
other transient radio signals. 



and the value of fh is in the range from to 1, depending 
on the configuration of the planets' collision. 

The distribution of the fragments by sizes for the de- 
struction of gravity-dominated bodies was presented in 
7], which with good accuracy has the power-law form 
n{D)dD cx D-^J^+^UD with the index l3 - 3.5-4. In this 
case, the index in ^ is 7 = (/3-H3)/3 ~ 2.17-2.33. So, 
the mass integral is dominated by the small masses, 
and the upper limit is not important. Then the number 
of fragments with the mass greater than m, is 



Nf{> m) 



(7 - 2)fhMh 



(24) 



Pieces of solid rock from the cores of the giant plan- 
ets could be formed only in collisions with large relative 
velocities, which is possible only in exceptional circum- 
stances (on the orbits very close to stars or after the 
supernova explosion). At smaller scales, where the role 
of gravity is low, the destruction should be accounted by 
the traction or granular material in a non-uniform shear 
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stress. The nature of fragmentation is determined by 
the strength of the substance and the spectrum of inho- 
mogeneities. Fragmentation due to shear stress occurs 
when the stress exceeds the breaking point of the sub- 
stance. In addition to the shear stress the shock waves 
are generated in the cohision, causing further cracking, 
and an additional characteristic scale is the wavelength 
of sound. The mass spectrum of the fragments may be 
associated with a spectrum of initial inhomogeneities in 
the core, if this spectrum was not lost due to the melting 
of the substance in the planet's interior. 

We now consider the thermal balance of the MH frag- 
ments on the orbit around a star after the matter re- 
lease in the collision of giant planets. The upper MH 
layers in Jupiter have the temperature T ^ 10* K. Im- 
mediately after the collision the ejected fragments will 
experience heating due to the transfer of the kinetic en- 
ergy into heat. Indeed, the kinetic energy of matter is 
Eq ~ If the ejected material would be heated 

uniformly, then its temperature would rise by the small 



amount AT ~ Eq/{I^MC) 



J{2C) - 4 X 10-3 K, 



where C is the specific heat capacity of MH. As C, we 
took the heat capacity of liquid hydrogen from [2^ 

C(T) = 6.86+0.66xl0-'*T+0.279xl0-^r2 kJ kg-^R-i 

(25) 

near T - 10"^ K. 

Suppose, the heating was completed and some frag- 
ment was ejected into space. It losses the heat by 
the thermal radiation and gets it from the Sun. Let 
the fragment has the spherical shape with the radius r, 
then its thermal energy is £' = {'^n / ?>)r^ pC {T)T , where 
/9 ~ 1 g cm"'^ is the MH density. This energy evolves 
according to the equation 



dE 
'dt 



(1 - A)LQ7rr^ 



fp{m,t), (26) 



where a2 is the Stefan-Boltzmann constant, Lq is the lu- 
minosity of the Sun, rorb is the radius of the orbit, A is the 
Bond albedo of the MH sample, and the unknown func- 
tion fp{m,t) describes the growth of the thermal energy 
due to the energy release during the phase transition of 
MH into the dielectric state. Then removing the pressure, 
after a period of metastability, the MH transformes from 
the metallic phase back into the dielectric, and the energy 
of the order of £ =290 MJ kg^^ releases [22|. It is possible 
that the transition goes from the surface into the sample's 
volume, or the transition happens in the whole volume 
at once or gradually, but in many local areas. In the case 
of a gradual transition the function fp{m,t) can be esti- 
mated as fp em/r, where t is the time of metastability. 
In the case of sharp transition fp ^ em5{t~ti — T), where 
S{t) is the Dirac delta function, and ti is the moment then 
pressure was removed. Remind that the theories of the 
MH metastable state were developed, but this question 
has not been clarified quantitatively either theoretically 
or experimentally, so the metastability must be regarded 
as the hypothesis with the free parameter — the time- 
scale of metastability. In view of these uncertainties in 



the estimates below we consider such the times t — U <^ r 
when the function fp{m,t) in the equation (j26p can be 
neglected. 

In the stationary case dE/dt = the temperature of 
the fragment is 



(1 - A)L, 



167r(T2r 



,2 

orb 



^fl-AV/\ r„.b X-V2 



82 



I- A 
"02" 



/ '-orb \ 

V5.2 AU/ 



K. 



(27) 

By the similar method, through the balance of absorption 
and emission of energy, the temperature of the asteroids 
is calculated usually. If the albedo of the MH is not 
extremely close to unity, the equilibrium temperature, as 
in the case of the asteroids between Jupiter and Mars, 
will be Te ^ 100 K. Note that in other systems, the 
luminosity of the stars may differ significantly from the 
luminosity of the Sun, and the characteristic temperature 
at the same distance from the star may be different. 

From (pS]) we obtain the cooling time from the initial 
temperature Ti to the certain temperature Tf 



At 



P7_ J{C 

3^72 J 



rpdC\ 

''^'dT. 



(28) 



The temperature of the fragment is only asymptotically 
approaching with At — >■ oo, so it's convenient in prac- 
tice to put the final temperature for example Tf = 1.2Te, 
then will serve as the definition of the cooling time. 
If heat capacity is given by then for Ti ^ 10* K, 

and for the interval Te ~ 10 — 100 K we have from the 
with a good approximation the expression 



At « 0.9 



P 



Ig cm 



V 1 cm/ 



VlOO K 



days. (29) 



For this characteristic time the sample is cooled to the 
temperature of 20% higher than the temperature (|27p . 
The mass cooled during the At is equal to 



An ^ 



5.7 



P 



1 g cm 



At 
1 day 



Te 



100 K 



(30) 

If the MH metastability time is greater than (|29|) . the 
cooling can move the fragment to the region of the phase 
diagram in which the MH is more stable. Otherwise, 
with less cooling time or the lack of stability regions, the 
fragment will pass into the gaseous phase and disperse 
in space, and this transition will be accompanied by the 
heating of the sample with the considerable release of 
energy e =290 MJ kg^^ 



VIII. RARE TYPES OF METEORITES 

With the destruction of the inner regions of giant plan- 
ets the two main types of meteorites could be formed: the 
destruction of the rocky core gives the stone fragments, 
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and the destruction of the MH areas under certain condi- 
tions may lead to the ejection of the MH fragments. As 
stated above, the question of MH metastability is still 
open, so the existence of the hydrogen meteorites should 
be considered as a hypothesis. 

Stony meteorites from the cores of the giant planets 
differ from the ordinary stony (and iron) meteorites in 
the sense that in the interior of giant planets, the sub- 
stance has been exposed to the much higher pressures 
than in the depths of the terrestrial planets or during 
the formation of small planetesimals, where the pressure 
was negligible. In addition, in the cores of the giant plan- 
ets the rock was in contact with dense hydrogen medium, 
which could also be reflected in the chemical composition 
of the meteorites. To clarify the possible chemical com- 
position and structure of the stony meteorites formed in 
the collisional destruction of the giant planets, more re- 
search is needed. 

Fragments of the MH can approach the Earth in 
the case of prolonged metastability at high temperature 
(about 100 K) or inside the meteorites, comets or aster- 
oids, been in the shade. In a worst case, when the MH 
has not the metastable state, the fragments of the heavy 
elements can still reach the Earth, which in a certain 
amount were mixed with MH in the giant planets and 
ejected during collisions. Perhaps such fragments can be 
identified in space according to their structural features 
and anomalous chemical composition, but the question 
about the characteristics of these fragments also requires 
investigation. 

Let us now consider the possibility of the MH mete- 
orite with a quite long metastability to be in the near 
space. By the ability to reflect the radio waves the MH 
fragments are not much different from the conventional 
iron meteorites. In optics, they have probably a differ- 
ent reflectance (albedo) than that of stone or ice bodies. 
However, for the direct optical observations the pieces of 
the size of the asteroid are necessary, which appears to 
be very rare. A more promising possibility is to identify 
the pieces of MH on their flight in the Earth atmosphere. 
During the combustion of hydrogen in atmospheric oxy- 
gen a trail of water vapour will be formed. 

It was shown in the Section |ll] that at the early stages 
of the planetary systems formation the collisions of the 
giant planet are possible. Suppose that at least one such 
encounter has been in the solar system with the mass 
ejections of MH fhMh - AMj/100 w IO-^/^Mq (see 
SectionlTVl , and the distribution of the fragments is given 
by We also assume that the fragments of MH re- 

mained metastable for about 5 billion years. Then the 
flow of fragments of MH onto Earth at the present time 
can be estimated as j{> m) ^ Nf{> m)vorh/V, where 
Vorh 30 km s~^ and V irL^a is the previously used 
characteristic volume of the protoplanetary disk. Numer- 



ically, for example for 7 = 2.25, we obtain 

.(>m)-l/.(^^j (ikij ^ ■ 

(31) 

For //i ^ 1 this value is much higher than the obser- 
vational constraints. Effects that could reduce the flow 
of debris on the ground are the following: the evapora- 
tion of the fragments, the collision of planets at larger 
distances from the Sun and the inability of the frag- 
ments to migrate to the Earth's orbit. But the sim- 
plest explanation is the small mass ejection in the col- 
lision //i ^ 1, or the absence of the collisions in the 
solar system. The large mass of MH in the solar system 
Mh ^ Mj/lOO « 10"^Mo probably contradicts the data 
on the movement of the planets and the limits of dark 
matter in the solar system. Indeed, the total mass of 
ordinary asteroids in the asteroid belt is estimated to be 
4% of the Moon mass. It is therefore necessary to as- 
sume //i <C 1, but even in this case, the flow ([511) can be 
significant. 

Now we estimate the flux of the MH fragments from 
other stars, swept into interstellar space by the tidal 
forces, as it was discussed in the Section |Vl In this case 
we have 2 x 10^ destructions with the same mass ejec- 
tions fhMh, but now distributed over the volume of the 
disk of the Galaxy Vd ^ irR^Hd, where Rd ~ 10 kpc 
Hd ^ 1 kpc, and the characteristic velocity dispersion of 
stars in the disk will enter instead of Worb- It should be 
noted that due to gravitational focusing the concentra- 
tion of fragments near the Sun will be somewhat higher 
than in the interstellar space. But since all the charac- 
teristic velocities are ~ 30 km s"""^, this effect gives the 
correction of the order of unity. For these numerical val- 
ues we have the coefficient 1.5 x 10~^^ yr""'^ km~^ in ([5T|) 
instead of the ~ 1 yr"-*^ m~^. I. e., in this case the frag- 
ments with mass > 1 kg fall every ^1000 years onto the 
entire surface of the Earth. 

Finally, we consider the intermediate case, when the 
destruction of the planets was near one of the neighbour- 
ing stars, or near the previous generation of stars in the 
vicinity of the solar system. Such MH fragments could 
be present now in the solar system. In this case, the 
mass fhMh is distributed over the average volume n~^, 
attributed to a single star. The corresponding coefficient 
in ([31]) is 3 X 10^^ yr^^ km~^, i. e. in the best case 
//j ~ 1 on average 1.5 fragments with mass greater than 
1 kg falls to the Earth per year. 

Thus, the flow of fragments of MH outside the solar 
system is very low. Detection of the fragments is possi- 
ble only when the collision occurred in the solar system, 
but the ejected mass of MH is much less then 10~^Mo, 
either there were barriers for the fragments to achieve 
the Earth's orbit. Nevertheless, the estimate ([51]) shows 
that there are good chances that the fragments may be 
present in an amount that is accessible for observations. 
Note that the hypothesis of the large MH fragment as 
a possible explanation for the Tunguska meteorite was 
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proposed by M. N. Tsynbal in 1983 and was developed 
by V. S. Mar kin. 

Of the particular interest are the small fragments of 
MH, which can reach the surface of the Earth by losing 
velocity (or entered the atmosphere with low velocity), 
but not burned completely. These fragments will evap- 
orate hydrogen from their surfaces, and its combustion 
can create volume glow with visual effect similar to the 
ball lightning. Ability to explain the mechanism of ball 
lightning by the burning of MH pieces was proposed in 
[27] . Ball lightning is actually the only macroscopic phe- 
nomenon, whose nature is not yet clear, see the review 
of observations and models in [28|, as well as separate 
works [13, [sO], [3l|. There are probably several differ- 
ent phenomena, which are called the ball lightning. For 
example, some of them may even be the processes in the 
visual cortex of the brain caused by electromagnetic pulse 
of the usual linear lightning fs^ . Flux of MH fragments 
can fall to Earth, producing the effect of the ball light- 
ning. Intensive surface combustion can create some lift 
due to the upward flow of hot gas. With regard to the 
MH fragments this issue requires further study. 

IX. CONCLUSION 

The collisions of giant planets in the process of the 
planetary system formation are fairly frequent, there are 
^ 0.1 — 1 collisions per each planetary system. In the col- 
lisions, the planets merge, loss the energy of the relative 
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